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ABSTRACT 

This article deals with a nonrelativistic quantum mechanical study of a dynamical sys- 
tem which generalizes the isotropic harmonic oscillator system in three dimensions. The 
Schrodinger equation for this generalized oscillator system is separable in spherical, cylin- 
drical, and spheroidal (prolate and oblate) coordinates. The quantum mechanical spec- 
trum of this system is worked out in some details. The problem of interbasis expansions of 
the wavefunctions is completely solved. The coefficients for the expansion of the cylindrical 
basis in terms of the spherical basis, and vice-versa, are found to be analytic continuations 
(to real values of their arguments) of Clebsch-Gordan coefficients for the group SU(2). 
The interbasis expansion coefficients for the prolate and oblate spheroidal bases in terms 
of the spherical or the cylindrical bases are shown to satisfy three-term recursion relations. 
Finally, a connection between the generalized oscillator system (projected on the z-line) 
and the Morse system (in one dimension) is discussed. 







Introduction 

The purpose of this paper is to study the quantum mechanical motion of a particle in 
the three-dimensional axially symmetric potential 

2 2 z z 2 x z + y z 

where fl, P, and Q are constants with Q > 0, P > —4, and Q > 0. In the last decade, 
this potential (including the case P = 0) has been the object of numerous studies [1-12]. 
The Schrodinger and Hamilton- Jacobi equations for this generalized oscillator potential 
are separable in spherical, cylindrical, and spheroidal (prolate and oblate) coordinates. 
In the case when P = we get the well-known ring-shape oscillator potential which was 
investigated in many articles [T], |3], f|, [?|, |8| in recent years as a companion of the Hartmann 
potential [4, 6, 8, 13-17]. If P = Q = we have the ordinary isotropic harmonic oscillator 
in three dimensions. 

The plan of this article is as follows. We solve the Schrodinger equation 

= EV, (2) 

with 

H=~(d axB +dyy + d tg )+V, (3) 

in spherical coordinates (in the second section) and in cylindrical coordinates (in the 
third section). (The constant U and the reduced mass are taken to be equal to 1. In 
the whole paper, we use \& to denote the total wavefunction whatsoever the coordinate 
system is; the wavefunctions \I/ in spherical, cylindrical, and spheroidal coordinates are 
then distinguished by the corresponding quantum numbers. Note also that we use s to 
denote the fraction | in the following.) In the fourth section, we determine the interbasis 
expansion coefficients between the cylindrical and the spherical bases. The fifth and sixth 
sections deal with the generalized oscillator system in spheroidal coordinates. In particular, 
the prolate and oblate spheroidal bases are expanded in terms of both the spherical basis 
and the cylindrical basis. Two appendices close this article. The first one is devoted to the 
bi-orthogonality of the radial wavefunctions (in spherical coordinates) for the generalized 
oscillator system. The second appendix concerns a connection between the Smorodinsky- 
Winternitz system (that is a basic component for the generalized oscillator system) and 
the Morse system. 

The generalized oscillator system constitutes a pending part to the generalized Kepler- 



Coulomb system studied in Ref. |TT|, [18| . The latter two nonrelativistic systems gener- 
alize two important paradigms in quantum mechanics, namely, the oscillator system and 

concern 



the Kepler-Coulomb system. The main results of this work and the one in Ref. []n_ 
the separability in spheroidal coodinates as well as the SU(2) approach and the three-term 
recursion relations for the interbasis expansion coefficients. 

The authors are very pleased to contribute to this memorial volume in honour of 
Jean-Louis Calais. Professor Jean-Louis Calais achieved, among other important works, 



an original job [19| on the derivation of the SU(2) Clebsch-Gordan coefficients by the 
(Lowdin) projection operator method. We are glad to present here a work where an 
analytic continuation of SU(2) Clebsch-Gordan coefficients plays an important role in the 
analysis of interbasis expansions. 
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The use of spheroidal coordinates is now well established in quantum chemistry [2(| . 
There exist now powerful techniques ]2T| for evaluating (angular and radial) prolate 
spheroidal wavefunctions from differential equations. It is hoped that this paper will 
shed some new light on expansions of spheroidal wavefunctions. 

Spherical Basis 

The Schrodinger equation (fj) in spherical coordinates (r,9,<p) for the potential ([!]), 



i.e. 



V 



n 2 



2 T + 2 r 2 cos 2 fl ' 2 r 2 sin 2 #' 



+ 



Q 



may be solved by seeking a wavefunction \l/ of the form 



V(r,6,<p) = R(r)Q(9) 



(4) 



with m G Z. This amounts to find the eigenfunctions of the set {H, L z , M} of commuting 
operators, where the constant of motion M reads 



M = L + 



P 



+ 



Q 



cos 2 9 sin 2 9 



(5) 



(L 2 is the square of the angular momentum and L z its ^-component). We are thus left 
with the system of coupled differential equations: 



1 



(M — A)Q 

A R 



-dJr 2 d r ) + 2E- tt 2 r 2 - 



0. 
0. 



(6) 
(7) 



where A is a (spherical) separation constant. 

Let us consider the angular equation (|G]). By putting 0(#) = /(#)/Vsin6>, we can 
rewrite Eq. (0) in the Poschl- Teller form: 



f 1 b 2 - i 

d ee + A + 

, 4 cos 2 9 



sin 2 9 



f = 0, 



''4 



Q 



In the case where b > s, the angular potential is repulsive for 9 — |. In this case, the 
9 domain is separated in two regions (9 G ]0, f[ and 9 G ]| , 7r[) and the "motion" takes 
place in one or another region. Furthermore, in this case Eq. (||) corresponds to a genuine 
Poschl- Teller potential. In the case where < b < s, we can call the angular potential 
an attractive Poschl- Teller potential. When b = s, i.e., P = 0, we get the well-known 
ring-shape oscillator potential |I|, |, f§ 0, §. The solution Q(9) = Q q (9; c, ±6) of Eq. (§) 



(for both < b < s and b > s), with the conditions 9(0) = @(f ) 
be (cf., H H) 

9(9) = N q (c, ±b){sm9) c {cos9) s±b P^ ±b) {cos29), 



0, is easily found to 
(9) 
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with q e N, where P&'P) denotes a Jacobi polynomial. Then, the constant A is quantized 

as 

A q (c, ±b) = (2q + c ± b + s)(2q + c ± b + 3s). (10) 
The normalization constant N q (c, ±6) in (|9|) is given (up to a phase factor) by 



@ q i@ q sm.9d9 = -S q > q - 



'IV, 



This leads to 



N q (c,±b) 



(2q + c±b + l)q\V(q + c ± b + 1) 
\ T(q + c + l)T(q±b+l) 



:i2) 



Note that only the positive sign in front of b has to be taken when b > s while both the 
positive and negative signs have to be considered for < b < s. 

Let us go to the radial equation (|7|). The introduction of (|T0|) into (0) yields an 
equation that is very reminiscent of the radial equation for the three-dimensional isotropic 
oscillator except that the orbital quantum number / is replaced by 2q + c ± b + s. The 
solution R(r) = R nrq (r; c, ±6) of the obtained equation, in terms of Laguerre polynomials 

R(r) = N nrq (c, ±b)(Vnr) 2q+c±b+ %- snr2 Ll q r +c±b+ \nr 2 ), (13) 
with n r G N. In Eq. (|T3p, the radial wavefunctions R nrq satisfy the orthogonality relation 

/ R n ' rq Rn rg r 2 dr = 5 <nr (14) 

" 

[cf., Eq. ([n])] so that the normalization factor N nrq (c, ±b) is 



N nrq (c, ±b) 



2Q 3s n r \ 



\ T(n r + 2q + c±b + 2)' 



(15) 



The normalized total wavefunction \l/(r, 9, if) = ^/ nrqm (r,9,ip;c,±b) is then given by 
Eqs. (||), (0), (|12D, (p!3|), and (|1^). The energies E corresponding to n r + q fixed are 



E n (c, ±b) = Q (2n + c ± b + 2), 



(16) 



with n = n r + q. Equation (|H]) shows that, for each quantum number n, we have two 
levels (for +b and —b) in the < b < s region and one level (for +6) in the b > s region. 
Note that this spectrum was obtained through a path integral approach in |2j, B] for the 



b > s case and in [|T^] for the general case (see also Refs. |jTD|, [TT| ) . 

In the < b < s region, for the limiting situation where b = s~, i.e., P = 0~, we have 
for the separation constant A: 



A q (c, +s) = (2q + c + l)(2q + c + 2), 



A q (c,-s) = {2q + c){2q + c+ 1). 



(17) 
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Then, by using the connecting formulas [24] 

(A) n +i 



^2n ( X ) 



xP( A "^(2x 2 -l), 



between the Jacobi polynomial P^ a '^' and the Gegenbauer polynomial C^, we have the 
following odd and even angular solutions (with respect to cos 9 i— ► — cos 9) 



Q q (9;c, +s 



(4g + 2c + 3)(2g+ 1) 



\| 27rr(2g + 2c + 2) 
x2T(c + s)(sin 0^(7^ (cos 0), 



6,(0; c 



(4g + 2c+l)(2g)! 



7 ^ 27rr(2g + 2c+ 1) 
x2T(c + s)(sin0) c C' 2 c + s (cos0). 



(19) 



Let us introduce (a new orbital quantum number) I and (a new principial quantum number) 
N through 



I — [ml 



iV — |m| 



2q + 1 for the + sign 
2q for the — sign J ' 



2n + 1 for the + sign 
2n for the — sign 



(20) 



Note that N = 2n r + l both for the + and — signs. Then, the separation constant [Eq. (|T 
and the energy [Eq. ([H])] can be expressed as 



A q (c, ±s)=A l (5) = (l + 5)(l + S + l), 
E n (c, ±s) = E N (S) = Q(N + 6 + 3s), 



(21) 



respectively, where 



Q + m 2 



m . 



Thus, the two parts of the energy spectrum for the signs ± correspond now to odd 
(for +) and even (for — ) values of N — \m\. In terms of N, I, and 5, the functions 
R nrq (r;c,±s) = R m (r;6) [cf, Eq. (|IJ)] and 6,(0; c, ±s) = 6 im (0;<5) [cf, Eqs. (0) and 



(|19|)1 can be rewritten as 



Rm(r; 5) 



2^(^)1 

2 

4 



(22) 



e, m (M) = 2 w+ *r(|m| + <5 + 



(21 + 25 +!)(/- Hj! 
AJ 27rr(/+ |m| + 25 + 1) 

x(sm^)l m l +<5 ^ l +f +s (cos0). (23) 

Equations and (p3| compare with the corresponding results for the ring-shape oscillator 
in H], [3|]. Note that (|23"D was given in terms of Legendre functions in Refs. [[J and || and 
was studied in details in Ref. ||. 

In the b > s region, for the limiting situation where b = s + , i.e., P = + , we have only 
odd solutions. In other words when P — > + , the eigenvalues and eigenfunctions of the 
generalized oscillator do not restrict to the eigenvalues and eigenfunctions, respectively, of 
the ring-shape oscillator. This fact may be explained in the following manner. To make 
P = in the wavefunction \l/„ rgm (r, 8, tp; c, +b) amounts to changing the Hamiltonian 
into a Hamiltonian corresponding to P = and to introducing an unpenetrable barrier. 
(Another way to describe this phenomenon is to say that for very small P, the potential 
V is infinite in the 9 = | plan and equal to the ring-shape potential only for P = 0.) This 
phenomenon is known as the Klauder phenomenon |25| . 

A further limit can be obtained in the case when 5 = 0, i.e., Q = 0. It is enough to 
use the connecting formula 

( _n\\m\ 

p^(x) = L -^ r r(M + s)(i - , 2 ) s Hc|* s (x) 



between the Gegenbauer polynomial and the Legendre polynomial P z . In fact for 
Q = 0, Eq. (|23D can be reduced to 



0,Je;O) = ^-l) , "' , 



!i±i(^)!pH (cos9) , 

\ 2 (l + \m\y. 

so that Qim(0; 0)e imip / V2tt coincides with the usual spherical harmonic Yi m (6, cp) (up to a 
phase factor, e.g., see |2l|). The wavefunctions g (6 l ; c, ±b)e imv / V2n may thus be consid- 
ered as a generalisation of the spherical harmonics. 

Cylindrical Basis 

In the cylindrical coordinates (p, <p, z), the potential V reads 



T . ft 2 , 2 2 , PI Q 1 

Equation with this potential, admits a solution \I/ of the form 



*{p,<p,z) = R{p)Z{z)-=, (24) 

V Z7T 

where m G Z. In other words, we look for the eigenfunctions of the set {H, L Z ,N} of 
commuting operators, where the constant of motion is 

N = D zz + ^, (25) 

z z 
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D zz being the zz component of 



D 



the so-called Demkov tensor [27] for the isotropic harmonic oscillator in R 3 . It is sufficient 
to solve the two coupled equations 



-d p ( P d p ) + 2E P - n 2 P 2 
p 



(N-2E Z )Z = 0, 
Q + m 2 ' 



R = 0, 



(26) 
(27) 



where the two cylindrical separation constants E p and E z obey E p + E z = E. The solutions 
(p, if, z; c, ±6) of d26)-23) lead to tne normalized wavefunction 



V 2 , «) = ^n p (p; c)Z p (z; ±6)- 



/ 2tt' 



(28) 



where 



^n p (p; c) 



and 



Z p (*;±6) = (-If J W e-^\Vnzy± b L±\Q; 



\p±b + 

with n p G N and p G N. The normalization of the wavefunction ( p8| ) is ensured by 



(29) 



/■oo 

/ Rn' p Rn p pdp = 8 n ' p n p , 
j 



ZpiZpdz — —dp/ p . 



Furthermore, the constants E p and E z in ( p26| - [27l ) become 



E p (n p ,c) = Q(2n p + c+V 



E g (p,±b) =n(2p±b + l). 



(30) 



Therefore, the quantized values of the energy E are given by (jig) where now the quantum 
number n is n = n p + p. As in the second section, the sign in front of b in Eqs. ([2~8])- (|3~0l) 
may be only positive when b > s. When < b < s, both the signs + and — are admissible. 
In the < b < s region, in the limiting case where b = s~ , due to the connecting 



formulas pi 



n 2n+1 (x) = (-l) n 2 2n+1 n\xL+ s (x 2 ) 
n 2n (x) = {-l) n 2 2n n\L~ s {x 2 
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between the odd H.2 n +i and even 7i 2 n Hermite polynomials and the Laguerre polynomials 
L^ s , we immediately have 



Z p (z;+s) = (-Y r -!— H 2p+1 (Vttz), 
W J2^+\2p+l)\ * 



Z p (z;-s) = (^j ^=== H 2p (Vttz). 



2 2 P(2p)\ 



Introducing (a new quantum number) n.3 such that = 2p + 1 for the + sign and = 2p 
for the — sign, we obtain 



The energy is then given by (|21|) where N = 2n p + n 3 + \m\. Note that the spectrum in 
the case b = s~, which corresponds to the ring-shape oscillator system, was obtained in 
Refs. [1-3]. 

In the b > s region, in the limiting situation where b = s + , we get only the odd solution 
of the ring-shape oscillator system. 

Connecting the Cylindrical and Spherical Bases 

According to first principles, any cylindrical wavefunction (pi]) corresponding to a given 
value of E can be developed in terms of the spherical wavef unctions (^) associated to the 
eigenvalue E (see also Ref. flllf )■ Thus, we have 

n 

^n pPm = T, W n P (c,±b)^ nrqm , (31) 
q=0 

where n p + p = n r + q = n. In Eq. Q3"T|), it is understood that the wavefunctions in the 
left- and right-hand sides are written in spherical coordinates (r, 9, if) owing to p = r sin 9 
and z = rcos9. The dependence on e imip can be eliminated in both sides of Eq. (|3lD . 
Furthermore, by using the formula L^(x) ~ (— l) n x n /n\, valid for x arbitrarily large, (|31|) 
yields an equation that depends only on the variable 9. Thus, by using the orthonormality 



relation (11), for the quantum numbers q, we can derive the following expression for the 



interbasis expansion coefficients 

W^c, ±b) = (-l)^^(c, ±6) E« np (c, ±6), (32) 

where 



BUc, ±b) 



(2q + c ± b + l)(n - q)\q\T(q + c±b+ l)T(n + q + c±b + 2) 
\| (n - p)\p\T(q + c+ l)T(q ±b + l)T(n -p + c+ l)T(p ±6 + 1) 

E q np (c, ±b) = 2 J\sm9) 2n - 2p+2c (cos6) 2p+1±2b P ? (c ' ±&) (cos 20) sin9d9. 



(33) 
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By making the change of variable x = cos 29 and by using the Rodrigues formula for the 
Jacobi polynomial f24[ 



-I) 



d r 



x) 



(1 _ x) -^i + x )-f>_ [(i _ x y+n {1 + x y + n^ 



Eqs. 



W«(c, ±6) 



2 n n\ 

lead to the integral expression 

(-l)P rl 



d q 



x, 



(2g + c ± 6 + l)(n - q)\T(q + c±b+ l)T(n + q + c±b + 2) 



(34) 



' ^ p!g!(n - p)\T(q + c + l)r(g ± 6 + l)r(n - p + c + l)T(p ±6 + 1) 
for the coefficient W% p (c, ±6). Equation ([34]) can be compared with the integral represen- 



tation 26 



(abaf3\cy) 



(2c+l)(J + l)!(J-2c)!(c + 7)! 



2 J+i 



^ (J - 2a)!(J - 26)! (a - a)!(a + a)!(6 - /?)!(& + /3)!(c - 7)! 

'-^(l + z)- 7 - 26 ]^ 
(a6a/3|c7) of the compact 



X^— / (l- x )*~«(l + x )>-f> [(l- x 



1 



(with J = a + b + c) for the Clebsch-Gordan coefficients C ( 
Lie group SU(2). This yields 



cy 
aa\i 



WUc, ±6) = (-I)""' (a 6 a/3|c , « + /?), 



rz ± 6 



a 



n + c c ± 6 

Co = gH — ^— , 



' u 2 

^ =F 6 
a = p — , p 



(35) 



n + c 



Since the quantum numbers in Eq. (|35| ) are not necessarily integers or half of odd integers, 
the coefficients for the expansion of the cylindrical basis in terms of the spherical basis 
may be considered as analytical continuation, for real values of their arguments, of the 
SU(2) Clebsch-Gordan coefficients. The inverse of Eq. fl31~D, namely 



* W = E^(C,±6)* 

rippm 

p=0 



(36) 



follows from the orthonormality property of the SU(2) Clebsch-Gordan coefficients. The 
expansion coefficients in (|36"D are thus 

W% q (c,±b)=W% p (c,±b). 

Note that in order to compute the coefficients W% (c, ±6) through (| 



]), we can use the 

3-^2(0, 6, c; d, e; 1) representation [|(| of the SU(2) Clebsch-Gordan coefficients. 

We close this section with some considerations concerning the limiting cases (P = 0, 
Q 7^ 0) and (P — 0, Q — 0). It is to be observed that the passage from (P 7^ 0, Q 7^ 0) 
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to (P = 0, Q 7^ 0) needs some caution. Indeed for b = s , Eq. ( j35| ) can be rewritten in 
terms of the quantum numbers N, I, and n 3 as 

W% p (c, ±s) = (-l)^(a 6 «/3|c ,a + /3), 

A" — \m\ — s ± s 



a 



4 



iV + Iml — s =F s 5 

= — ■ h - 

4 + 2' 

_ 21 - 1 5 

2n 3 — N + \m\ — sis 

4 ' 

-2n 3 + AT + Iml + s ± s 5 
= — + -. 

H 4 2 



By using the ordinary symmetry property [26 



(aba0\cy) = (~l) a+b - c {ab, -a, -0\c, -7) 

and the Regge symmetry [p6|| 

, . . /a + 6 + 7 a + 6 — 7 a — 6 + a — a — b — a + 0. ^ 
(aba0\c^) = I , , , \c,a-b 



(37) 



in Eq. ([37]) with the sign + and by using the ordinary symmetry property 

{aba0\c-f) = (-l) a+b - c (ba0a\ C1 ) 
in Eq. ([T?]) with the sign — , we get 

W% p (c, ±s) = W l Nmn3 {5) = (a o b o a0\c o , a + /?), 

N+\m\ 5 1 iV-|m|-l 21 - 1 5 

a o = ~ A + 2' b o = 4 ' c o = + 2' 

AT + |m| - 2n 3 5 n 2n 3 — A" + |m| — 1 
« = : 1" 7:, P 



4 2' ' 4 

As a conclusion, when 6 = s~ we have an expansion of the type 

^Nmn 3 {p, <P, Z;$) = Y1 W NmnM ^Nlm(r, 6, (p; 5), (39) 
I 

where the summation on I goes, by steps of 2, from |m| or |m| + 1 to iV according to 
whether as N — |m| is even or odd (because N — I is always even). Equations (|38|)- (|39l) 
were obtained in Ref. H for the ring-shape oscillator system. Finally, the case P = Q = 
can be easily deduced from (|38l)-(j39"l) by taking 5 = 0: we thus recover the result obtained 
in Refs. [ |2l| P9"| for the isotropic harmonic oscillator in three dimensions. Note that in the 
case P = Q = 0, the expansion coefficients in Eq. ( |3~9"D become Clebsch-Gordan coefficients 
for the noncompact Lie group SU(1,1) (cf., Ref. fl30|| ). 
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Prolate and Oblate Spheroidal Bases 

SEPARATION IN PROLATE SPHEROIDAL COORDINATES 

The prolate spheroidal coordinates (£, rj, (p) are defined via 



R 



2 
R 



(e-l)(l-^)sin^, 



(with 1 < £ < oo, — 1 < r\ < 1, and < (p < 2tc), where R is the interfocus distance. As is 
well-known |[31|| , in the limits where i? — > and R — > oo, the prolate spheroidal coordinates 
reduce to the spherical coordinates and the cylindrical coordinates, respectively. In prolate 
spheroidal coordinates, the potential V reads 



p 



+ 



Q 



ev 2 (e-m-v 2 ) 

By looking for a solution ^ of Eq. (0), with the potential ([40]), in the form 

pinup 



(40) 



2vr 

with m G Z, we obtain the two ordinary differential equations 

P 



(41) 



d^e - ik - 



Q + m 2 ER 2 n Vt 2 R 4 2 2 



+ 



-r- 



16 



1 + 



, c2 



d,,(l -T] 2 )d ri 



Q + m 2 £i? 2 
1 - ;/-' 2~ '' 16 // 



^ = +A(i?)Vi, 



fi 2 i? 4 



'02 



-A(i2)V: 



2- 



(42) 



(43) 



where X(R) is a separation constant in prolate spheroidal coordinates. The combination 
of Eqs. ( f42"D and fl43|), leads to the operator 



A 



1 



v %(e-m+ed v (i-v 2 )dn 



+ 



+ 



£ 2 — rf 

e+v 2 -i (n fl s 

(£2 -1)(1- ^2) W °WJ 

e 2 + ^ 2 



n 2 R A 



erf + P 



£2^2 



after eliminating the energy E. The eigenvalues of the operator A are X(R) while its 
eigenfunctions are given by (f4i~D . The significance of the (self-adjoint) operator A is to be 
found in the connecting formula 



d2 

A = M + —N, 



(44) 
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where M and N are the constants of motion (j^) and (p5|) . The operator A is of pivotal 
importance for the derivation of the interbasis expansion coefficients from the spherical 
basis or the cylindrical basis to the prolate spheroidal basis. In particular, it allows us 
to derive the latter coefficients without knowing the wavefunctions in prolate spheroidal 
basis. (In this respect, credit should be put on the work by Coulson and Joseph [j3^] 
who considered an operator similar to A for the hydrogen atom.) Therefore, we shall 
not derive the prolate spheroidal wavefunctions ipi and ip2 which could be obtained by 
solving Eqs. ( f42[) and ([43l) . It is more economical to proceed in the following way that 
presents the advantage of giving, at the same time, the global wavefunction \l/(£, 77, <£>) = 
Vi ¥1 R-i c i an d the interbasis expansion coefficients. 



INTERBASIS EXPANSIONS FOR 

THE PROLATE SPHEROIDAL WAVEFUNCTIONS 

The three constants of motion M, N, and A, which occur in Eq. 
satisfy the following eigenequations 



Mm 



n r qm 



Uppm 



A q {c, ±6)* 
2E z (p, ±b)m nppm , 



, can be seen to 

(45) 
(46) 



and 



A$ nkm = X k {R)m nkm (47) 

for the spherical, cylindrical, and prolate spheroidal bases, respectively. [In Eq. ([47]), the 
index k labels the eigenvalues of the operator A and varies in the range < k < n] The 
spherical, cylindrical, and prolate spheroidal bases are indeed eigenbases for the three sets 
of commuting operators {H, L z , M}, {H, L z , N}, and {H, L Z ,A}, respectively. We are now 
in a position to deal with the interbasis expansions 



Vnkm = J2 U nk(R;C,±b)y npPm 
p=0 
n 

^nkm = J2 T nk(R;C,±b)^ nrgm 
q=0 



(48) 
(49) 



for the prolate spheroidal basis in terms of the cylindrical and spherical bases. 

First, we consider Eq. (|48"D. Let the operator A act on both sides of (f48"|). Then, by 
using Eqs. ([44"1), (f4~6D , and (|4~7|) along with the orthonormality property of the cylindrical 
basis, we find that 



where 



\ k (R)-^E z (p, ±b) 



U p nk (R;c } ±b) 



J2u£ k (R;c,±b)M$, 

p'=0 



M {± ) 
pp' 



00 r2ir poo 



JO 



r(±) 



pdpdipdz. 



(50) 



(51) 



The calculation of the matrix element Mp pl can be done by expanding the cylindrical 
wavefunctions in (|5T|) in terms of spherical wavefunctions [see Eq. (31)] and by making 
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use of the eigenvalue equation for M [see Eq. (f45|) 1. This leads to 



i n 

KP = o E ±b)WZ p (c, ±b)WZ p ,(c, ±6). (52) 



To calculate the sum in Eq. fl52|) , we need some recursion relation for the coefficient 
W% p (c, ±6) involving p — 1, p, and p + 1. Owing to Eq. (|35|), this amounts to use the 
following recursion relations [j33fl : 



[-a(a + 1) - 6(6 + 1) + c(c + 1) - 2a/3](abaP\cj) 

= yj(a + a) (a - a + 1)(6 - /3)(b + (3 + 1) 
x(a, 6, a — 1, /5 + 1)07) 

+ ^(a - a)(a + a + 1)(6 + f3){b - (3 + 1) 

x(a,6,a + l,/5 - l|c7). (53) 

Then, by introducing Eq. (^) into Eq. (|52|) and by using the ort honor mality condition 

E(a6a/3|c7)(a6a'/5'|c 7 ) = S^S^, (54) 



c,7 



we find that M pp ) is given by 



M$ = 2\p(p ±b)(n-p + l){n + c-p + l)] s V, P -i 
+ [s (c =f 6 + s) (c =f 6 + 3s) + 2(p + \){n - p) + 2(p ± 6)(n + c - p + l)]8 p > v 

+2[{p + l)(p+l±b){n-p){n + c- p)] s 5 p , !P+1 . (55) 

Now by introducing (|55"D into (0), we get the following three-term recursion relation 

[(p + 1) {n — p) + {p ± 6) (n + c — p + 1) 

+i(c T 6 + a)(c T 6 + 3s) + ^E z ( Pl ±6) - \\ k (R)]U p nk 

+ [(p+l)(p+l±b)(n-p)(n + c- p)] s Kt l 
+ [p(p±b)(n-p + l)(n + c-p + l)] s U p nk 1 = (56) 

for the expansion coefficients U^ k = U^ k (R; c, ±6). The recursion relation (^) provides us 
with a system of n + 1 linear homogeneous equations which can be solved by taking into 
account the normalization condition 

The eigenvalues Xk{R) of the operator A then follow from the vanishing of the determinant 
for the latter system. 

Second, let us concentrate on the expansion ([19|) of the prolate spheroidal basis in 
terms of the spherical basis. By employing a technique similar to the one used for deriving 
Eq. 0, we get 

[\ k {R)-A q {c,±b)]T« nk {R-c,±b) 

= ^Y.Ti(R;c,±b)N$, (57) 

Z q'=0 
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where 



qq' 



JO JO 



mi rqm Nm nrq , m r 2 sin edrdOd^p. 



The matrix elements N qq , can be calculated in the same way as M^ v , except that we must 



use the relation [26 



pp 



x • 



, nh R \ \ _[ c 2 (2c+l)(2c-l) 

1 P| ° 7j " [{c 2 - 7 2 )(-a + b + c)(a - 6 + c)(a + 6 - c + l)(a + 6 + c + 1) 

(c - 7 - l)(c + 7 - l)(-a + 6 + c-l)(a-6 + c-l)(a + 6- c + 2)(a + 6 + c) 

(c- l) 2 (2c-3)(2c- 1) 

. , n . , (a — /?)c(c — 1) — 7a(a + 1) + 76(6 + 1) . , „. 

x (a6a/3 c - 2, 7) - - — , \ - — -(aba(3\c - 1, 7 

c(c-l) 



and the orthonormality condition 



^2(aba/3\c r y)(aba{3\c''y') = <5 C / C 5y 7 , 

a,/3 

instead of Eqs. (|53| ) and (|54|) . This produces the matrix element 



AT (±) p , 1 ,\ 2g(g + l) + (c±6)(2g±6+l) ^ 
qq nK ' ' (2q + c±b)(2q + c±b + 2) qq 
-2fiK +1 (c, ±b)8 ql>q+l + A q n (c, ±&)V,,-iL 



(58) 



where 



A q n (c,±b) 



q(n — q + l)(q + c ± b)(q ± b)(q + c)(n + q + c ± 6 + 1) 
(2g + c ± 6) 2 (2g + c ± 6 - l)(2g + c ± 6 + 1) 



Finally, the introduction of (58) into (p7|) leads to the three-term recursion relation 



p2 

A fc (i*) - A q (c, ±b)-—E n (c, ±b) 



+QR Z 



2q(q + 1) + (c±6)(2g±6 + 1) 
X (2q + c±b)(2q + c±b + 2) 

A q n + \c,±b)T q t l + A q n (c,±b)T q ^ 



rpq 

1 nk 







(59) 



for the expansion coefficients T% k = T^ k (R; c, ±6). This relation can be iterated by taking 
into account the normalization condition 



T,\TMc,±b)\ 2 = l. 

q=0 

Here again, the eigenvalues Xk{R) may be obtained from the vanishing of the determinant 
of a system of n + 1 linear homogeneous equations. 
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LIMITING CASES 

Putting b = s~ , i.e., P = 0~, in the matrix element 
we have 



with Q ^ and by using 



qq' 



2A l {5)-2{\m\+5) 2 -l 
En {6) (2/ + 25-l)(2/ + 25 + 3)^ 
-2Q[A l + 2 (5)5 V)l+2 + A l N (tyi>j-2), 



where 



/_(/_ - + 25)(l+ + 26- 1)(N - / + 2)(N + 1 + 25+1) 



4(2/ + 25- l) 2 (2l + 25- 3)(2Z + 25+1) 
(with l± — I ± |m|) and finally we get the following three-term recursion relation 



wm 4^ t „ 2A l {5)-2{\m\+5f-l 
X k (R) - A t (5) - Y E ^ 6 \2l + 25-l)(2l + 25 + 3) 



T l Nk (R;5) 



+nR 2 [A l + 2 (5)T l N \ 2 (R; 5) + A l N (5)T l Nk 2 (R; 5) 

for T l Nk (R;5) = T l Nk (R;c,±s). By analogy it is easy to obtain a three-term recursion 
relation for the interbasis expansion coefficients U^ 3 k (R] 5) = U^ 3 k (R; c, ±s). We get 

[(2n 3 + l)(iV - n 3 + 5 + 1) + (\m\ + 5) 2 - 1 + ^-(2n 3 + 1) - A A (i?)]C/^(it!; 5) 

+ [(n 3 + l)(n 3 + 2)(JV - \m\ - n 3 )(N + \m\ - n 3 + 25)} s U n N 3 k +2 {R- 5) 
+ [n 3 (n 3 -1)(N- \m\ -n 3 + 2)(N + \m\ - n 3 + 25 + 2)} S U^ 2 (R; 5) = 0. 

Consequently, when b = s~ we have the expansions [cf., Eqs. fl48|) and (|49|)1 

N 

n 3 
N 

Nlmi 



for the ring-shape oscillator. The summations on I and n 3 go, by steps of 2, from \m\ or 
\m\ + 1 to N and from or 1 to AT — |m| according to whether as iV — \m\ is even or odd 
(because iV — I and N — |m| — n 3 are always even). 

The next limiting case 5 = 0, i.e., Q = 0, is trivial and the corresponding results for 
the isotropic harmonic oscillator agree with the ones obtained in Ref. f3Hfl . 

Finally, it should be noted that the following two limits 



lim U* k (R;c,±b) 
lim T« k (R;c,±b) 



W% k (c,±b), 
W q nk {c,+b) 



furnish a useful means for checking the calculations presented in the fourth and fifth 
sections. 
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SEPARATION AND INTERBASIS EXPANSIONS FOR 
THE OBLATE SPHEROIDAL WAVEFUNCTIONS 

The oblate spheroidal coordinates (£,rj,ip) are defined by 

R 
2 



x = ^-\f (£ 2 + - r] 2 ) cosip, 



R r=o, 

.'/ = ■g-V (f + -^ 2 )sin^, 

z = 

(with < £ < oo, — 1 < 77 < 1, and < ip < 2ir), where R is the interfocus distance in the 
oblate spheroidal coordinate system. As in the prolate system, in the limits R — > and 
R — > 00, the oblate spheroidal coordinates give the spherical and cylindrical coordinates, 
respectively |31], j34f . 

The potential V, the Schrddinger equation, the oblate spheroidal constant of motion 
A, and the interbasis expansion coefficients for the oblate spheroidal coordinates can be 
obtained from the corresponding expressions for the prolate spheroidal coordinates by 
means of the trick: £ — ► — z£ and R — > %R. 

Spheroidal Corrections for the Spherical and Cylindrical Bases 

As we have already mentioned, the spheroidal system of coordinates is one of the 
most general one-parameter systems of coordinates which contains spherical and cylindri- 
cal coordinates as some limiting cases. Accordingly, the prolate spheroidal basis of the 
generalized oscillator as R — * and R — > 00 degenerates into the spherical and cylindrical 
bases that can be treated as zeroth order approximations in some perturbation series. The 
three-term recursion relations for the expansion coefficients of the prolate spheroidal basis 
in the cylindrical and spherical bases, which have been obtained in the fifth section, may 
serve as a basis for constructing an algebraic perturbation theory, respectively, at large 
(R 3> 1) and small (R ^ 1) values of the interfocus distance R. Thus it is possible to 
derive prolate spheroidal corrections for the spherical and cylindrical bases. 

THE CASE R < 1 

Let us rewrite the three-term recursion relation (|59D in the following form 



[X k (R) - A q (c, ±b) - nR 2 B« n (c, ±b)]Tl k 
+{IR 2 [Al +1 (c, ±b)T q n t 1 + Al(c, ±b)T^ 1 ] = 0, (60) 

where 

Bl{c,±b) = l -(2n + c±b + 2) 

2q(q + l) + (c±b)(2q±b+l) 
X (2q + c±b)(2q + c±b + 2) 

The zeroth order approximation for the separation constant Xk(R) and the coefficients 
T^ k (R;c,±b) can immediately be derived from the recursion relation (|60|). Indeed, from 
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Eq. (pOf) , we obtain 



\im \ k (R) = A k (c, ±b), 

R — >0 

l im n T nfc(^; c, ±6) = 5 kq , 
it — >u 

so that, for the wavefunction, we have 

lim ^nfcm(^, 7), tp; R, c, ±6) = ^„ rfcm (r, 0, ip, z; c, ±b). 

R — >0 

As is seen from these limiting relations, the quantum number k, labeling the spheroidal 
separation constant and being (according to the oscillation theorem ||31| ) the number of 



zeros of the prolate angular spheroidal function ^(v) m t ne interval — 1 < r\ < 1, turns 
into a spherical quantum number determining the number of zeros of the angular function 
(|9p. It is clear that this fact is a consequence of the independence of the number of zeros 
of the wavefunction on R. 

In order to calculate higher order corrections, we represent the interbasis coefficients 
T% k (R;c,±b) and the spheroidal separation constant X k (R) as expansions in powers of 
ttR 2 : 

oo 

T« k (R; c, ±6) = 5 kq + £ ({lR 2 y, (61) 

3=1 

oo 

\ k (R) = A k (c, ±b) + £ \® ({lR 2 y. (62) 

3=1 

Substituting Eqs. (^l|) and ( |62]) into the three-term recursion relation flBOD and equating 
the coefficients with the same power of R, we arrive at the equation for the coefficients 
Zj£? and \ k j) 

A(k-q)(k + q + c±b + l)T k i i ) = 
-A^(c, ±b)T^l + Bi(c, ±b)T^ 

-A« n (c, ±b)T<?-V - £ Xt^S- (63) 
t=o 

Equation (|63"D with the initial condition T kq = 8 kq and the condition T^> = 5jo arising 
in the standard perturbation theory [|35| allow us to derive a formula expressing Ajjf for 
j > 1 through the coefficients T^ 1 ^ and T kk ±i- 

A? = ~A k n +1 (c, ±b)Tj?-V + S*(c, ±b)T^ - A k n (c, ±b)T^\. (64) 

This gives a possibility to determine, step by step, the coefficients X k ^ and in Eqs. (|6T|) 
and (0). As an example, let us write down the first and second order corrections in (|62| ) 
for X k (R) and the first order correction in (|6lD for T^ k (R; c, ±6). It follows from Eq. (|64|) 
that 

Al 1} =^(c,±6), 
= -^ +1 (c, ±6)T« +1 - A k n (c, ±6)^ 
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and Eq. (^3|) for j = 1 results in 



T (i) = 



A k n (c,±b) 



A k +\c,±b) 



k « ~ A(2k + c±b) 6q ' k - 1 + A(2k + c±b + 2) 6q ' k+1 - (65) 

Thus, for the spheroidal separation constant, with an accuracy up to the term (QR 2 ) 2 , we 
get 

X k (R) = A k (c, ±b) + nR 2 B*(c, ±b) 

Q 2 R 4 \ A k n {c, ±6) 2 _ A k n +1 (c, ±bf 
2k + c±b ~ 2k + c±b + 2 



Introducing ( |65|) into (|6l"l) and then using (49) for the expansion of the prolate spheroidal 
basis over the spherical one, we get the following approximate formula 

^nfcm(£, 77, (p\ R, c, ±6) = ^ nkm (r, 6, tp; c, ±6) 
Q 2 R 4 \ A h n {c,±b) 



Aj + \c,±b) 
2k + c±b + 2 



THE CASE R > 1 

Now let us consider the case R 3> 1- The three-term recursion relation (|56|) can be 
written as 



where 



[D p n (c, ±b) + ^E z (p, ±b) - ^}U p nk 
+ [Cr 1 (c, ±b)Ul + k l + C*(c, ±b)Ul~ l ] = 0, 



C£(c,±6) = [p(p±b)(n-p + l)(n + c-p+l)] s , 

D p n {c,±b) = {p+l)(n-p) 
+ {p ± b){n + c — p + 1) 



+ -(c^6 + s)(c^6 + 3s). 



It follows from Eq. 



that 



(66) 



l im MS = ij5;(]fc,±6), 

R^oo 

For i? 3> 1, the interbasis expansion coefficients and the spheroidal separation constant 
are developed in negative powers of QR 2 : 



U p nk {R- c, ±6) = S kp + £ tfg>(ni*V, 



(67) 
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Substituting Eqs. ( |67]) and fl68[) into Eq. (|66"D , we get 

+C£(c, ±b)U^l\ - i^A?^^ = o. 
4 t=l 

Using the conditions llj^ = 5k p and Ujjj) = 5jo, one easily obtain 



(68) 



0+1/ 



rtf-l) 



(i-i) 



r(i-l) 



(69) 



(70) 



Equations ( p9[ ) and ( |70|) completely solve the problem of determining the expansion coef- 
S°andl7£> 



ficients A! and . For instance, we have the approximate formulae 
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^(c,±6) 2 -Q fc+1 (c,±6) 2 



(fii? 2 ) 2 

+ C n( C > ±fe )^,fc-l,m(P ) ^Z;C,±&) 

-C^ +1 (c, ±6)^ nifc+1 , m (p, z; c, ±6) 
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Appendix: Bi-Orthogonality of the Radial Wavefunctions 

Besides the orthonormality relation (|nj) in the quantum numbers n r for the function 



R 



nrq , we also have an orthogonality relation in the quantum numbers q, viz., 



x(±) 

J qq' 



RnwRnrqdr = 2q + c±b + l 5q ' q ' 



(71) 



for a given value n' r + q' = n r + q of the principal quantum number n. The proof of (|7T[) 
is as follows. In the integral in Eq. (ITT]), we replace the two radial wavefunctions by their 



expressions (|lq). Then, with the help of the formula pq] 

POO 

/ Q- cx x a - x Llicx)LHcx)dx = 
Jo 

(7 + l) m (A-a + l) n r(q) 
m\n\c a 

X3-F 2 (— m, a, a — A; 7 + 1, a — A — n; 1), 
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we arrive at 



7 (±) 



n 



T(q' + q + c±b + l) 
T(2q' + c±b + 2)T(q-q' + l) 



X, 



(n- q')\T(n + q' + c ± 6 + 2) 
(n-g)!r(n + g + c±6 + 2) 

x 2 F 1 (-q + q', q + q' + c± b + 1; 2g' + c ± 6 + 2; 1) 



By using the Gauss summation formula p4[ 

r(c)r(c-a-6) 



2 Fi(a, b; c; 1) 



we can rewrite (|72|) as 



7 (±) 



r(c-a)r(c-6) 



ft 

+ + c±6 + l 



(72) 



(n- q , )\T(n + q' + c±b + 2) 
\ (n- q)\T(n + q + c±b + 2) 

xlTiq-q' + l^iq'-q + l)]- 1 . 

This completes the proof of Eq. ( ffT] ) since [T(q — q' + l)T(q' — q + = 8, 

Appendix: The Smorodinsky-Winternitz and Morse Systems 

The Morse system with the potential 

V M = V (e- 2ax - 2e~ ax ) 

can be connected to the dynamical system with the potential 

fl 2 2 PI 
Vsw = Y z + --. 

(The latter potential may be considered as a one-dimensional component of the so-called 



Smorodinsky-Winternitz |i7[ |38], |39| potential. The potential Vsw was investigated by 
Calogero PD|.) 



The Schrodinger equation for the Morse potential Vm, i-e., 

~d xx + 2E- 2Vo{e~ 2ax - 2e~ ax )] $ = 



(73) 



admits a discrete spectrum (with E < 0) and a continuous spectrum. For the discrete 
spectrum, by making the change of variable 



-sy 



y = ax, 1/6R, z = e 
and the change of function 



z G R 
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in Eq. 



), we get 



f = o, 



(74) 



where 



A 



Equation (|74]) has the same form as Eq. (^H) for z > with 

E 7 = 4A 2 , = 2A, P = 



8E 

2 



(7 



Therefore, we must consider two admissible regions for the energy E: (i) -32E > a 2 and 
(ii) < -32E < a 2 . 

In the case (i), by employing the energy formula (|30|) for E z , we obtain that E is 
determined by the relation 



V^2E 



A-(p + s), p = 0, [A 



In Eq. (|75|), [a;] stands for the integral value of x. As a result, we have 



E = -Vn 



1 , ^ 

1 - ^ (p + s 



p = 0,l, 



[A -5] 



(75) 



(76) 



Equation (|76|) is in agreement with the well-known result according to which the discrete 
spectrum of the Morse system has a finite number (here [A — s] + 1) of energy levels with 
the condition A > s. 

In the case (ii), we have 



V^2E 



±[A-(p + s)], 



which has no solution for p G N. 

The connection just described between the Morse and Smorodinsky-Winternitz systems 
can be used also to deduce the wavefunctions of one system from the wavefunctions of the 
other. For instance, from Eq. (p9|), we immediately get the normalized solution ip(x) = 
^;A)offD: 



V>p(*;A) = (-l)*(2A) 



A-p 



apl 



r(2A-p) 



xe 



-Az 2 ^2A-2p-l£2A-2p-l 



with 



[ (2Xz 2 



p = 0, 1, • • • , [A - s] 



(77) 



Our result (77) differs from the one of Nieto and Simmons |4]J (by the fact that the factor 
p! in fl77|) is 2p- A in Ref. H). 
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